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I Abstract 

j^jq' In this paper we study the future asymptotics of spatially homogeneous Bianchi type II cosmologies 

with a tilted perfect fluid with a linear equation of state. By means of Hamiltonian methods we first 
find a monotone function for a special tilted case, which subsequently allows us to construct a new 
set of monotone functions for the general tilted type II cosmologies. In the context of a new partially 
gauge invariant dynamical system, this then leads to a proof for a theorem that for the first time gives 
a complete description of the future asymptotic states of the general tilted Bianchi type II models. 
\ The generality of our arguments suggests how one can produce monotone functions that are useful 

■ for determining the asymptotics of other tilted perfect fluid cosmologies, as well as for other sources. 

^ ' 1 Introduction 

Spatially homogeneous anisotropic perfect fluid models have during the last decades been successfully 
studied using a dynamical systems approach. The book [1] summarizes most of the presently known 
^ , results about the so-called non-tilted perfect fluid cosmologies, while the more general 'tilted' perfect 

\^ ' fluid models have been primarily investigated more recently [2] — 114) . 

In all of the papers investigating tilted models, the analysis has rested on standard techniques from 
dynamical systems theory. Most of the results concern the identification of fixed points and a subsequent 
linear stability analysis of these points. In order to get a grip on the global aspects of the solutions, an 
effective tool is the use of monotone functions. However, such monotone functions are hard to find, and 
in most of the previous works on tilted models the monotone functions were obtained by brute force, trial 
and error, and luck. It would therefore be desirable to have a more systematic method to seek and find 
such monotone functions. 

For non-tilted spatially homogeneous perfect fluid models, virtually all known results crucially rely on 
the existence of conserved quantities and monotone functions. These have turned out to be connected to 
the existence of certain symmetries, intimately associated with conservation laws such as the preservation 
of the number of particles in a fluid element, and the so-called scale- automorphism group [15]. Although 
not necessary, the symmetries and associated structures were, to a large extent, found by means of 
Hamiltonian techniques, see ch. 10 in [T] and [T5]. One aim of this paper is to illustrate that one can 
fruitfully use similar methods that previously have been applied to non-tilted models to tilted ones. To 
do so, we will consider an example — the tilted Bianchi type II models. 
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2 DYNAMICAL SYSTEMS DESCRIPTION OF TILTED BIANCHI TYPE II COSMOLOGIES 2 

The tilted perfect fluid Bianchi type II models have been analyzed before as a dynamical system in [5] . 
In that paper, as well as in the present one, the perfect fluid was assumed to obey a linear equation of 
state characterized by p/p = w = const, where p and p are the pressure and energy density with respect 
to the rest frame of the fluid, respectively; special cases of interest are dust, w = 0, radiation, w — ^, 
a stiff fluid, w = 1, while a cosmological constant A can formally be regarded as a perfect fluid with 
w = —1, however, in this paper we will consider the range — 1 < it; < 1. 

In [5] the fixed points of the system were found and their linear stability properties were studied. It 
was observed that the future stability of the fixed points depended on the equation of state parameter w, 
but that a future stable fixed point existed for all w in the range — 1 < w; < 1 under consideration. On 
the basis of the linear analysis and numerical computations it was conjectured that the future linearly 
stable fixed points were global attractors in the full state space. This conjecture was corroborated by 
means of monotone functions for some ranges of w. However, the complete picture was not fully drawn, 
mainly because of the lack of a set of sufficiently restrictive monotone functions. In this paper we will 
use methods that previously have been applied to the non-tilted models and find a set of new monotone 
functions that are sufficiently restrictive to determine the future asymptotics of all tilted solutions, thus 
filling in the missing gaps in the conclusions drawn in [5]. 

The outline of the paper is as follows. In Section [2] we give a new partially gauge invariant dynamical 
system for the general tilted Bianchi type II models, derived in Appendix|Xl Then, based on the structure 
of monotone functions obtained for more special models by means of Hamiltonian methods, presented 
in Appendix [B] we obtain the new monotone functions in Section [3l which leads to a new theorem that 
describes the future asymptotic states of all tilted Bianchi type II cosmologies. Finally, in Section |4] we 
conclude with a discussion about the structure of the monotone functions, why they exist, and why one 
can hope to expect similar structures in other models. 

2 Dynamical systems description of tilted Bianchi type II cos- 
mologies 

In [16j the so-called conformally Hubble-normalized orthonormal frame equations are given in full gener- 
ality. These are in turn specialized to the spatially homogeneous Bianchi case in Appendix [X] and then 
to the presently studied Bianchi type II models with a general tilted perfect fluid with a linear equation 
of state. Then a subsequent set of new variables, invariant under frame rotations in the 23-plane (see 
Appendix |Al , yield the following state vector and dynamical system. 

State vector. 

S-(S+,E,E2^s2^r!k,«'), (1) 
where we treat E^, E^, and as variables, but where we have refrained from giving them new names. 

Evolution equations: 



Constraint equation: 





= -(2 - q)E+ - 3E2 + 41]k + ^(1 + w)G^^v^n , 


(2a) 


E' 


= -(2 - g)E - 2V3E^ + Vst^ + ^(1 + w)G+\^n , 


(2b) 




= -2(2 - q - 2\/3E)E^ 


(2c) 


(E^)' 


= -2[2 - g - 3E+ + %/3E]S2 , 


(2d) 




= 2(g - 4E+)r!k , 


(2e) 




= 2GI^ (1 - v^)[3w - 1 - E+ - V3t]v^ . 


(2f) 




/(S) = 4S2nk - (1 + wfG^^v^n^ = . 


(2g) 



In the above equations is given by the Gauss constraint 

= 1 - E^ - r^k , (3) 
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where 

= + f]2 + £2 ^ (4) 

The condition Q > in combination with ^ yields < + 'S'^ + & + t'^ + fl]^ < 1. The deceleration 
parameter q is given by 

q^2T,^ + ^G^^[l + 3w+{l-w)v'^]n, (5) 

while 

Gi = 1 ± ww^ ; (6) 

finally, ' denotes differentiation with respect to a dimensionless time parameter r, determined by dr = 
Hdt, where t is the clock time along the congruence normal to the spatially homogeneous hypersurfaces. 

We now give a brief description of the invariant subsets and fixed points of the system ([2]), which 
is analogous to the analysis given by Hewitt et al. [5]. Note that although our system is invariant 
under frame rotations in the 23-plane it is not invariant under all rotations. Hence there exist multiple 
representations of solutions, for further comments on this, see [5]. 

Name Restrictions Dimension Interior/Boundary 



i) 


Non-tilted non- vacuum Bianchi type H 


= 


S2 


= 


4 


Boundary 


ii) 


Non-tilted non- vacuum Bianchi type I 


= 




= 


4 


Boundary 


iii) 


Vacuum Bianchi type H 




: n : 


= 


4 


Boundary 


iv) 


Vacuum Bianchi type I (Kasner) 




: n - 


= 


4 


Boundary 


v) 


Extreme tilt 


= 


1 




4 


Boundary 


vi) 


Orthogonally transitive Bianchi type H 


S2 = 


- 




4 


Interior 



Table 1: Invariant sets of the state space. The last column indicates if a subset is part of the boundary 
of the general tilted Bianchi type II models or if it is an interior subset. 

Fixed points: 

i) The flat Friedmann solution, F: —l<w<l, 

S I = E = = ^2 = = t;2 = . 



ii) The Collins-Stewart solution, CS: — ^ < w < 1 



-(3i«-|-l), S = = ^2 = ^2 = 0, f7k = — (3w-l-l)(l-w). 



64 



iii) Hewitt's solution, H: | < w < 1, 



S^. = -(3^ + l), S = — (7^-3), E =0, E ^ 16(17.. -1) 

3(1 -u;)(5w + l)(3w- 1) ^ (3w-l)(7w-3) 

"k = TT-Z -T , V = 



4(17u>-l) ' {llw + 1){5W + 1) ' 

iv) Hewitt et al.'s 1-parameter set of solutions^ HL: w = |, < 6 = const < 1, 

f:^', E2.±6, E2,4(46+l)(8-35)^ 
+ 3 ' 3V3 27 ' 513 

^ _ (2b+ l)(17-86) ^2 _ 3(4b-H)(26-f 1) 
~ 171 ' " - (i7_86)(8-36) ■ 



^There is a misprint in [5] where the square root on the b in the Si expression of the line of fixed points, HL, has 
disappeared. 
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v) Hewitt et al.'s extreme tilted point, Het: —1 < w < 1, 

^1^1 ~o 4 100 _ 3 



^ 3 ' 3V^ ' 27 ' 513 ' 19 

The system also admits the following fixed point sets: the Kasner circle K*-^, for which = 0, the Kasner 
lines KL±, for which — const, and the extremely tilted Kasner circle Ket°, for which = \. These 
subsets reside on the Bianchi type I vacuum boundary, i.e., = 1, with E = S = 0, see [5], however, 
since these fixed points do not play a prominent role in this paper we refrain from giving them explicitly. 

Remark. Below we will refer to the relevant fixed point values for S^. and S by S+o and Eqi respectively. 

3 Future asymptotes in tilted Bianchi type II cosmology 

In what follows certain monotone functions will play a crucial role. Based on our results in Appendix \K\ 
andlHl we hence begin by deriving them. 

3.1 Monotone functions 

There are several auxiliary equations that are useful in the context of monotonic functions, see Ap- 
pendix |X| 

Auxiliary equations: 

n' = [2q - (1 + 3w) + (1 + w)(3w - 1 - E+ - V3t,)G^\^]n , (7a) 
Q' = -[2(l-g) + E+ + V3S]Q, (7b) 
= [2g- (l + 3w)]^', (7c) 

where 

Q^{l + w)G+-lvn, ^T-'-^-'''^ G+^n. (8) 

Since 

2q-{l + 3w) = 4I]2 - (1 + 3w){l - n) + {1 - 3w){l + w)G+\!'^n > , if - 1 < w < -1/3 , (9) 

as follows from (|47|) . ^ is a monotonically increasing function when —l<w< —1/3; henceforth we 
denote ^' in this interval of the equation of state by Afp- 
Before we continue, let us introduce some notation: 

= 1 - E+oS+ - SoS , (10a) 
ip^, = [S+o(S+ - S+o) + So(S - to)f , (10b) 
= [Eo(S+ - S+o) - S+o(S - to)f , (10c) 

where the subscript * henceforth denotes a specific fixed point, while S+o Si,nd Sq are the associated fixed 
point values for and S, respectively. In the following it is important that > 0, which can be seen 
as follows: 

0, - i [1 - S^o - + 1 _ _ E2 + (S+ - S+o)2 + (S - tof] > i [1 - E% - > 0, (11) 

where we have used the Gauss constraint ([3]) and 17 > 0, ilk > 0. 
For non-tilted perfect fluid models 



n'^n^-"^ 3(i-w)s+o „ 1 



Mcs:v-o = 0CS f^'"" = -VqS+)^ ' ™ " 8(1 - £%) ' = 8^^ + ^""^ ' ^° " ° ' ^^^^ 

is a monotonically increasing function. However, it is of interest to generalize this function by replacing 
with which is equal to f2 in the non-tilted case, i.e., 

Mcs = 0cif^" (13) 
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which leads to the following time derivative in the present fully tilted state space: 



cs 



(l-w) 



'PCS 



(1 ^+o)^+o 



^+0 



1 



^ + o (3 - 7w)(2S'^ + (1 + w)G+ v^n) 



, (14) 



and hence Mcs is monotonically increasing when —1/3 < w < 3/7. 

In Appendix |B] we derive a monotonic function for the tilted orthogonally transitive case 
which can be written as 







-(3+13«i) 



'(S") — ri^"^ • (15) 
It sometimes turns out to be the case that a monotone function for a given state space is also monotone 
in a more general state space in which the original is embedded in, at least for a limited range of the 
equation of state parameter, see e.g. [lOl [15]. We hence compute the time derivative for Mr in the full 
tilted case; this gives us 



(lnMH)'-0H' 



8 + |(91w — 31)(7w — 3) 4 



(16) 



and hence A/h is monotonically increasing when 3/7 < w < 5/9. 
The above monotonic functions all have the form 

i), and cki, a2, as, a4, /3 > 0. For the individual cases we have 

1^2 = as = , a4 = 1 , 13 = 2, (18) 
\2 — , as = m , a4 = 1 — m, /3 = 2 , (19) 
) , ^2 = 5(7u; - 3) , as = 3w - 1 , a4 = 4 , /3 = 3 + 13w . (20) 

Let us assume the form p!7|) in order to find a monotonic function for the range 5/9 < ui < 1. We 
obtain 

MHet=<^Hct(S^)^(S^)^X, (21) 

and hence ai = 4 , a2 = 10 , as = 9 , a4 = , /3 = 46, which leads to 



+ a2 H 


- as 


Mf : 








Mh : 





^Hot 



243^Het + 207^Hct + ^3(9^ 5) ^„i^^ _ ^^^^ 



(22) 



(luA/Hct)' 

and thus A/net is monotonically increasing when 5/9 < -u; < 1. 
3.2 Future asymptotic limits 

Let us denote the invariant set for the general tilted Bianchi type case, for which (1— tj^) v"^ J7 fik 7^ 0, 

by ^Gon, while we denote the orthogonally transitive case, for which (1 — v'^) i7 fit 7^ 0, by 5ot- 

A local analysis, as in [5], reveals that the state space S have fixed points as local sinks according to 
Table [21 which leads to the following bifurcation diagram: 



=-1/3 



> CS 



=5/9.6=1 



> Het. 



(23) 

It was conjectured in |5j that the above local results hold globally modulo a set of solutions of measure 
zero. Using the above monotone functions one can show that this is indeed the case, and also identify all 
the exceptional solutions. In the following theorem we show that all solutions that belong to Sqcu end 
up asymptotically at the above sinks: 



Theorem 3.1. For all x e Scen 



a;(x) 



F -l<u;<-l/3 

CS -l/3<w<3/7 

H 3/7<w<5/9 

HL w = 5/9 

Het 5/9 < w < 1 



4 DISCUSSION 



6 



Range of w 


Sink 


-1 < w < -1/3 


F 


-1/3 <w<3/7 


CS 


2>/7 <w< 5/9 


H 


w = 5/9 


HL 


5/9 < u; < 1 


Het 



Table 2: Sinks for iSgo 



The proof makes use of the Monotonicity Principle [1] I17| , which gives information about the global 
asymptotic behavior of solutions of a dynamical system. It is stated as follows: Let 0t- be a flow on R" 
with X being an invariant set. Furthermore, let M be a C"'^ function M : X — >■ R. Then if M is increasing 
on orbits, then for all x £ X 

uj{x) C {s e X\X I lim M{y) inf M} . (24) 

Proof. We make use that iScon is a relatively compact set and hence that every orbit in iScon has an 
w-limit point in Sqcu, moreover, the oj-limit set of every orbit in iSccn must be an invariant set. In every 
case inf = 0, and hence we only have to investigate the set where (InM*)' — 0. It turns out that in 

all cases the invariant set associated with (In Af*)' = is precisely the pertinent fixed point(s), which thus 
is the cj-limit of every orbit in Sqcu- Hence the proof, and the situations, is virtually identical to that of 
the non-tilted Bianchi type II perfect fiuid case given in [5] on p. 151. For the cases —l<w< —1/3, 
-1/3 < w < 3/7, 3/7 < w < 5/9 and w = 5/9, 5/9 < w < 1, one uses Mp, Mqs, Mh, and A/not, 
respectively. □ 

Remark. It follows that F, CS, H, Het attract a 4-parameter set of solutions, as does the line HL, however, 
in this case it follows from the reduction theorem, see e.g. [18], that each point on the line attracts a 
3-parameter set since the line is transversally hyperbolic. 

Corollary 3.2. For all x e Sqt 

{F -l<w<-l/3 
CS -1/3 < w < 3/7 . 
H 3/7<w<l 

Proof. This follows immediately from the previous proof, in combination with noticing the form for 
(IuA/h)' when E2 ^ 0. □ 

We hence have established that the local bifurcation diagram 



F ^^^^ CS -^^^ H (25) 



from [5] reflects the global features of the solution space of Sqt • 



4 Discussion 

In this paper we have found that there exists a collection of monotonically increasing functions that 
completely determine and describe the future asymptotics of tilted Bianchi type II models. Furthermore, 
they all take the form 

Af, =07^(S2)"i(E2)a.f^a3(^)a4^ (26) 

where /3 = 2(ai + a2 + 03 + 04), and ai, a2, 013, 0:4, /? > 0. Moreover, the time derivatives of the 
monotone functions McSi -^h, -^Hct, with nonzero S^q + Eg (and hence non-zero + <f*) all take the 
form 

(InAf,)' = (l)-^{A^(p., +Inv,), (27) 
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where A^ and are constants, and Inv* a function that vanishes on one of the invariant sets in Table 
[TJ Outside these invariant sets the function Inv* is positive only in a limited range of values for w. In 
the Mcs case Invcs is zero for the non-tilted subset for which = E'^ = 0; in the Mh case Invn is 
zero for the orthogonally transitive subset for which = 0; finally, in the Muct case Invuct is zero for 
the extreme tilt subset for which = 1. Clearly it is easier to find these monotone functions on these 
subsets first and then extend them to the larger state space, as we did in order to find the key monotone 
function Mh in Appendix [B] moreover, we also first computed A/not for the extreme tilt subset — 1, 
although we did not use Hamiltonian methods in this caseH 

The importance of the hierarchical structure of Bianchi cosmology, where we have systems with 
boundaries on boundaries, have been emphasized before, see e.g. [HJEniET] and references therein. Here 
we see yet another context for this observation, which suggests that one should first try to find monotone 
functions for subsets and then attack the case one is really interested in. Hence one should first identify 
subsets for a given state space and write them on the form Za — and then, if there exists a locally future 
stable fixed point on a given subset that admits subsets Za — attempt to find monotone functions of 
the form 

M,=<l>-f^l[z:-. (28) 

a 

Is there a deeper reason for why monotone functions like this should exist? The analysis of the 
non-tilted case in [15] suggests that the existence of these monotone functions are related to the scale- 
automorphism group. In the tilted case this group can be viewed as consisting of an off-diagonal special 
automorphism group and a diagonal scale-automorphism group. The off-diagonal special automorphism 
give rise to conserved momenta, if the underlying symmetry is not broken by source terms, and hence also 
to monotone functions, in a similar way as for the non-tilted models, see |T^. We have not discussed such 
monotone functions here since we did not need them for the future asymptotics, however, they could be of 
help for the much more difficult past asymptotic behavior. However, the off-diagonal automorphisms also 
have other dynamical consequences. It is because of the off-diagonal automorphisms we only had diagonal 
shear degrees of freedom in 0* in the present case, and it is because of this one would expect the present 
analysis to also be of relevance for other tilted models, and for other sources. Moreover, in the present 
case we have encountered a hierarchy of source subsets that is typical. It is the increasing complexity 
of the source that breaks the vacuum symmetry group associated with the scale-automorphism group, 
creating a hierarchy of monotone functions associated with different source subsets. This is also what 
one encounters in the non-tilted case jl5) . however, with an increasingly complex source this phenomenon 
seems to be even more pronounced! Hence a systematic attempt on the tilted models, or other sources, 
strongly suggests a deeper investigation into the dynamical consequences for the scale-automorphism 
group for the various relevant subsets. This is quite ambitious task and we have therefore refrained 
from doing it here; instead we have made use of the structures one can expect to arise from such an 
analysis without deriving all the details that completely determines all the monotone functions from the 
scale-automorphism group (for a hint of how this can be done, see the complete analysis of the class A 
non-tilted models from a Hamiltonian perspective in |15) ) in order to see if this is likely to be a fruitful 
project. The answer seems to be yes. 
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A Derivation of the dynamical system 

In spatially homogeneous cosmology the space-time is foliated by a geodesically parallel family of spatially 
homogeneous slices with a timelike unit normal vector see, e.g., [l][22l[23] and references therein. 

^If one wants to use Hamiltonian methods to deal with extreme tilt, one first has to observe that these models have the 
same equations as those associated with a source that takes the form of a perfect fluid with a radiation equation of state 
w = 1/3, but with a null vector field replacing the timelike 4- velocity. 
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A.l Perfect fluids 

Splitting the total stress-energy tensor Tab with respect to n" yields: 

Tab = pnanb + 2q(^a "b) +pKb + T^ab , (29a) 
p = n''n''Tab, qa = -ha''n^Tbc, p = ^ h'^' Tab , l^ab ^ h^a" hb^"" T,d , (29b) 

where hab = naTib + gab and A/^ab) = ha'^hb'^Acd — ■^hath'^'^Acd', p,p is the total energy density and total 
effective pressure, respectively, measured in the rest space of n". In this paper we consider a perfect fluid, 
which yields the stress-energy tensor: 

T'^''^ + (30) 

where p and p are the energy density and pressure, respectively, in the rest frame of the fluid, while -u" 
is its 4-velocity; throughout we assume that p > 0. Making a 3-1-1 split with respect to n", leads to 

=T{n'' + V''); Uav" = , T ^ {1 - v^)-'^/^ , (31) 

where is the three-velocity of the fluid, also known as the tilt vector; this gives 

p = T-^G+^p, q'' = {l + w)G+^ pv'\ p = wp+^{l-3w)qav\ T^ab = q{aVb) , (32) 

where G± — l±w , w — p/ p. 

A. 2 Orthonormal frame equations 

In Bianchi cosmology the metric can be written as 

''g = -N^{x°) dx^ ® dx° + 5-3-(a;°) {J + N' dx°) (x) {J + dx^) = 1, 2, 3) , (33) 

where {u;*} is a left-invarianl[l co-frame on G dual to a left-invariant spatial frame {ej}. This frame is a 
basis of the Lie algebra with structure constants C**-^, i-e., 

[ej, e.] - C^j. = (ej.,.„ n™^ + 2a[j S-/) e-^ or, equivalently, dJ - -i^.^ J A a;^ , (34) 

where we have decomposed the structure constants C'-^ as [23]: 

cSk = ^]k,r^ + ^^4^ (35) 

where a-- — ^^ej^-^ejpft n*'^n'"", — aio^ = \h[{n^i)'^ — n*-- n-' ■], where his a, constant group invariant 
that is unaffected by frame choices. The Bianchi models are divided into two main classes: The class A 
models for which = 0, and the class B models for which ^ 0. 

For simplicity we will set the shift vector to zero, i.e., N"^ — 0. This leads to the orthonormal frame 
''g = -{i^° f + (u;i)2 + (u;2)2 + (a;3)2^ and a frame {e„} that is dual to {w"}, a = 1, 2, 3, 

eo=iV-iA, e^=ej{x'')e. = ej{x'')e-;-^^ where ^"/'ej(x°) e;3^(x") = 5^^^(a;0) , (36) 

where g*-' (x") is the left-invariant contravariant spatial metric associated with gj-- {x^) (ef = ef{x^ , , x^)). 

Since the unit normal to the spatial symmetry surfaces n = Cq by definition is hypersurface forming, 
and since it is the tangent to a geodesic congruence due to spatial homogeneity, we obtain 

[ eo, ] C^oa ep ^ fj ep = -[H 5 J + aj + ej^ n'i]ep, (37a) 
[ea,e^] = C^a^e^ = [2a[„ (5^]'^ -f Cq/j^ '^*'''] , (37b) 

^See [l] p. 38 on the meaning of group invariant frames and their relation to orthonormal frames. 
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where H is the Hubble variable; aap is the shear associated with n; fl" is the Fermi rotation which 
describes how the spatial triad rotates with respect to a gyroscopically fixed so-called Fermi frame0 The 
relations (|36|) and (|37l) yield 



^=3 



(e^^) , n'^^ = .g-5e". n^^' , a„ = , (38b) 

where g is the determinant of the spatial metric gjj, and e"| is the inverse of Bq', i.e., 

g = det(ff53.) = (det(e"j))2 = (det(e„'^))-^ e^j e^^^ = . (39) 

The matter conservation equation "VaT""^ = for a perfect fluid with a linear equation of state yield 

{\npy = [1 + w)G^^[-3H + f^pv'^v'^ + 2aaV°'], (40a) 
v = GZ^il~v^) [3wH + fafsc^c^ -2wafsc'^ v]v, (40b) 



Cr 



- [Sj - c^A [// - v{ap + e^^5 n^^ c^)] , (40c) 



where ap + e^^^ rt'^'' c''' — C" c^; c''' , and where instead of the 3- velocity Va we have found it convenient 
to introduce v — \/vaV" > and the unit vector Cq = Va/v as variables. To obtain a more compact 
notation, we also introduced / = (/)' = Cq/ = N~^df /dx^ = df /dt, where t is the clock time associated 
with the normal congruence of the spatial symmetry surfaces (i.e. N = 1 for this parameterization). 

It is of interest to consider the particle number density h and the chemical potential jl, which, for a 
linear equation of state, can be defined as (see [231 [53] and references therein) 

p = /i = (1 + w) n'" . (41) 

Defining 

l = ngiT, (42) 

yields the evolution equation (In I)' = 2aav" = 2(aQc")w, and hence I is a constant of the motion 
whenever a^v'^ = 0, e.g. for the class A perfect fluid models. Another quantity of interest is Taub's 
spatial circulation 1-form |271 128) ta — p-Uai whose spatial components can be written as = jlT = 
fiTvCa, with the norm jiTv, which satisfies (In flT v)' = fapc^c^, which, together with (|40cl) yields 
ia — fiTv {fafi — V C^afi c-y) ■ This in turn leads to that t-- — e" ■ ta, where Cq (e"j) = —fp" e^i, obeys 
the equation 

i. = -{fiT)~^C%tf^l) . (43) 
A. 3 The Hubble-normalized dynamical systems approach 

In the conformal Hubble normalized approach one factors out the Hubble variable H by means of a 
conformal transformation which yields dimensionless quantities [291 116j . In the spatially homogeneous 
case this amounts to the following: 

{Y.ap,R'',N''f',Aa) = -|(a„;3,f^",n"'',aa), (f7, P, g„, H„^) = ^(p,p, ^„^) , (44) 

where we have chosen to normalize the stress-energy quantities with 3iJ^ rather than in order to 
conform with the usual definition of il; in the perfect fiuid case this leads to that 

= (l + iy)G;%17c", Q = {l + w)G+'^ vn, P = wn + \{l -■iw)vQ , H„^ = t; g C(„C/3) . (45) 

In addition to this we choose a new dimensionless time variable r by means of the lapse choice 
N — H~^. Since H is the only variable with dimension, its evolution equation decouples from the 
remaining equations for dimensional reasons: 

H' = -{l + q)H- q = 2Y? + + = (46) 



*The sign in the definition of is tfie same as in | 25II16| . but opposite of that in [l]. 
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where a prime denotes d/dr and where q is the deceleration parameter, obtained by means of one of 
Einstein's equations Gat = Tab — the Raychaudhuri equation (we use units c = 1 and SnG = 1, where c 
is the speed of hght and G is Newton's gravitational constant); in the perfect fluid case we obtain 

g = 2I]2 + i(17 + 3P) = 2I]2 + iG+^[l + 3u; + (1 - w)v^]n . (47) 

The remaining Einstein field equations together with the Jacobi identities yield the following set of 
equations, which we divide into evolution equations and constraints; 

Evolution equations 

E^^ = —(2 — q)^ai3 — 2e''''(Q S^)^ Rs — '^{a^) + 3X1^^ , (48a) 

(A^"'^)' = (3<7(5^(" - 2F^(")iV'')'^ , (48b) 

A'^ = Fj Ap , (48c) 

n' = (2q - 1) f7 - 3P + 2A„ Q" - E^^jH"/^ , (48d) 

v' = -GZ^ (1 - v^) [1 - 3w + E„0 c°'c^ + 2w {Ap c^^) v] v , (48e) 

c'a = [SJ - c^cP] [Fp^ - v{A0 + N^^ c^)] . (48f) 

Constraint equations 

= 1 - E^ - - , (49a) 

= (3(5,^^^ + 6,*'' 7V5^)E^^-3Q„, (49b) 

= Ap N^^ , (49c) 

where E^ = iEa^E"'^, = A^A", Fj = q5j - - ej^ i?^, and 

=%0> -2e^'(„iV^)^A5; n^^^B"„ + A^- B^p ^ 2N^^ N^^ ^ - N\ N^p . (50) 

The Gauss constraint (j49a| and the Codazzi constraint (j49b| will figure prominently throughout. 
It is of interest to note that 

g' - -[2 - <7 - F^p c" + 2(A„ c") v] Q . (51) 

Another quantity of interest, intimately connected with I, is defined via fiT = g~^/'^l] raising the r.h.s. 
with 1 + w and normalizing with H yields the quantitjH 

^ ^p-(i-t«)^-i^^ (ln4')'==2q-(l + 3u;) + 2(l + w)(^„c")?;. (52) 
A. 4 Bianchi type II 

For the Bianchi type II models we have A^ — 0, and in addition we can choose a spatial frame to be 
an eigenframe of the matrix N^p, with iVn ^ 0, while otherwise N^p = 0. This leads to that eq. (|48bl) 
yields i?2 = E31, i?3 = — E12, and that the Codazzi constraint (j49b[) implies 

wi = = ci . (53) 

Inserting the conditions of eq. ([55]) into eq. (|48f[) gives the following relation: 

= E12 C2 + E31 C3 = -i?3 C2 + i?2 C3 ^ CAB R^ = Q , (54) 

where A, _B = 2,3 and ^ab is the two-dimensional permutation that has £23 = 1 (hence cac^ — 1)- If 
follows from cab R^ = Q that Ra (x ca o< Qa, where the last relation holds when ^ 0; hence E12 
and E31 are linearly dependent and can be replaced by a single variable. 



^This quantity has appeared before in the hterature, e.g. in [5], where F ™) been denoted by , 
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We have the freedom to rotate in the 23-plane, which is expressed in the field equations as the 
freedom to choose i?i . To obtain a set of variables that are invariant under such rotations we introduce 
the following new shear variables 

E+^iS^^^-iEii, S = ^(Eab-S+'5ab)(c^c^-^^^''), (55a) 

E - ^(Eab - S+<5^s)e^c(c^c^ - ^S^'^) , = ^ (S?^ + ^li) ■ (55b) 



Hewitt et al. |5j make use of the freedom to rotate in the 23-plane to set C2 = 0, which yields that 
C3 = 1 and i?2 = = S31. This leads to a correspondence between the variables Si, E3 in ^ to our 
variables, when setting C2 — 0, according to S = — S_, = E3, = E^. 

Because of the existence of discrete symmetries one can simplify the analysis, e.g. the eigenvalue 
analysis, by introducing the following state vector S = (E+, E, E^, E^, fik, t'^), where fit = -^ii/12, and 
where can be obtained in terms of S via the Gauss constraint (|49ap . This leads to the dynamical 
system given in Section [21 

B Hamiltonian considerations and derivation of monotone func- 
tions 

B.l Hamiltonian considerations 

The scalar Hamiltonian is given by 

f( = 2Ng^n'' n\Gab - Tab) = 2iVgn" n^Gab - Tab) = N^K = iV(r + Ug + Ui) , (56) 

where we have defined N = Ng'^l"^, and where (T, C/g, Uf)=Qg H'^{-1 + E^, f^j,, n). By means of (|32l) . 
(|iT|l . and dH]), this leads to that 

Ui = 2gp ^ 21^+"" c,(i-"')/2 T^-^G+ , (57) 

while T and C/g depends on the model and the metric representation. 

In [30^ it was shown that the tilted orthogonally transitive Bianchi models exhibit a so-called timelike 
homothetic Jacobi symmetry. It was later realized that such symmetries are related to the existence 
of monotonic functions (Uggla in ch. 10 in jl], and [15]). Unfortunately the analysis in [3(T is quite 
cumbersome and we will therefore make a new derivation of the 'homothetic structure' and from this 
derive a monotone function. To do so, we need to connect the spatially homogeneous frame in (I33p 
with the orthonormal frame. This is done in two steps: (i) diagonalization by means of the off-diagonal 
special automorphism group, (ii) normalization by means of diagonal scaling. We hence write in the 
transformation (p6|) . i.e., — eQ*(x°)ej, as Bq* = (£'^^)q-' (5^^)j\ or e"-- — D"jS^i, where S^i, since 

it is assumed to be a special automorphism transformation, leaves oj and n'^ unaffected. In addition we 
define the new metric variables /3^,/3^,/3'^ via the matrix so that 

/exp(-/3i) \ 

iD-')o.' = exp(-/32) , (58) 

V exp(-;93)/ 

where l3" = /3"(a;0); hence 5^/2 = exp{f3^ + fi'^ + P^). 

To obtain the Hamiltonian for the tilted orthogonally transitive Bianchi models we choose a spatially 
homogeneous frame so that the line element can be written as (1331) with all structure constants being zero 
except n^^ = fii. In the orthogonally transitive case we can specify the spatially homogeneous frame so 
that g~-j in ([55)1 have one off-diagonal component, (7^2 j ^^'^ so that the perfect fluid velocity has a single 
non-zero component, Vf,. Hence we follow 130] and writeEl 

xO) 0\ 

(59) 




^There is a typographical error in eq. (2.61) in |30| ; the exponent should be —1 and not —1/2 of (m^^^) in the expression 
for S3. 
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It follows from ([aSb]) that 



n 



11 



exp(^i - - P^)ni . (60) 



Due to that 63 is associated with the off-diagonal special automorphism group, the associated momentum, 
which is proportional to (T12, is conserved [30]. Furthermore, note that (|43| yields the constant of the 
motion ^3 = ^:3= const. This is consistent with the Codazzi constraint, also known as the Hamiltonian 
momentum constraint, which linearly relates these two constants to each other. The constants fii, Z, ^3 
allows us to write T + Ug + Ui in ((55)) so that 

J{ = T + J7g + L/f = Td + J7c + C/g + f/f , (61) 

where Uc is the so-called centrifugal potential, which is proportional to the term in E^, where, 
furthermore, Uc, Ug, and U{ are all expressible in terms of /3", and no other time dependent quantities. 
Expressing NT^ in terms of the or the associated momenta tTq, , see [H] and also [31] , by means of (p8| 
and the transformations in this subsection, leads to 

NTd - 2iV-i g^5/3^/3* = jN g-'^n^ns , (62) 

where Q-ys is known as the minisuperspace metric for the diagonal degrees of freedom, which, together 
with its inverse 0°^^ is given by 





= -1 -1 ; = i -1 1 -1 . (63) 



The centrifugal potential Uc = 2g I]^2 = 2 5 cri2 re-expressed via the Codazzi constraint in terms 
of ^3 yields 

Uc oc exp[-2(/?i - /32)] . (64) 
By means oiUg — Qg ilk and (|60|) we find that 

[/g ocexp(4^i). (65) 

We can write = — 1 as = 1 — + jl^^ t", and by defining 

F = (1 + w)-^ t„ e = [l + U.)-2 g^ ^ , (66) 

we obtain 

= 1-F2-|-FF2'", (67) 

which allows one to, in general, implicitly express F^ in terms of F, i.e., F^ = F^(F). In the present case 
we find that F = {\ + w)-"^ l''^"' 5'" 5^3 ^. ^, ^ which yields 

F = (1 + w)-^ r^"" exp[w{l3^ + 13^) - (1 - w)(3^] . (68) 

It follows that 

iJf cxexp[(l-w)(/3i -f /32+^3)]0. (69) 
where 4> := F^~"'G'+ is a function of the particular combination w{/3-^ + 0^) — (1 — w)f3^ only. 



B.2 Derivation of monotone functions 

Based on the Hamiltonian for the diagonal degrees of freedom for the orthogonally transitive type II 
case we in this subsection derive a monotonic function that is of key importance for understanding 
the dynamics. In ch. 10 in [T| and in [15] it is shown that monotone functions are associated with 
'homothetic' symmetries of the potential, i.e., we require that there exists a vector c = c"9^q such that 
cU — c{Uc + Ug + Ui) = rU, where r is a constant. For this to be possible we require (i) that c0 = 0, 
and hence that 

c{wil3^ + 13^) - (1 - w)l3'^) = wic^ + c^) - (1 - wy = , (70) 
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and (ii) cUc ~ r Uc, cUg — r Ug, and cUf — r Ug, which due to the condition (i) yields cexp[(l — 
^2 _^ ^3)] ^ ^ gxp[(l - w){P^ + This leads to 



{1 - w){c^ + + c^) = r, 



(71) 



which yields = 3c^, — 4c^u;/(l — w); w.l.g we can choose — 1 ~ w, which gives 

(ci,C2,C3) = (l-u;,3(l-u;),4u;), r = 4(1 - m;) . (72) 
The causal character of this vector with respect to the metric Qafi is crucial |15l ; we obtain 

g^sc^ ^ -2{1 - w){3 + 13w) , (73) 

which is timelike if —3/13 < w < I. The above properties of the potential U implies that the model 
satisfies the criteria given in |15] for admitting a monotonic function given bj0 



M cx (c" 7r„)^exp 



(74) 



We now have to express M in the dynamical systems variables S . We do so in two steps by first using 
Sq,q — Sa, and then in the second step we go over to the presently used dynamical systems variables 
via the current gauge fixing (see the previous discussion about the correspondence between our variables 
and those used by [5j). Let us define 



then 



which yields 



and hence 



exp 



{Vg, Vt) = exp 



g-isc< /?' 



{Ug,Ui), 



exp 



J-C — exp 



Gjsc^ 



^Gjsc^ 
Vg + V, 



T 



oc 



T + V. + Vf^O. 



M cx 



(1 



(75) 
(76) 
(77) 
(78) 



since Vg/V{ = Q^/fl and 1 - T,'^ = fli, + Q. 

It follows from our definitions that — ^ttq (2 — S^^), ttq = 711+712 + ^3, see also [151 HH]- In 
combination with ([7^ we find that this yields 

^ (X 1 - i(l - «;)Ei - ^(1 - w)E2 - Iwi:^ = 1 - ^{7w - 3)S - i(l + 3^i;)E+ . (79) 
ttq o z o o 

Next we need to solve for Vf in terms of the state space variables (E+, E, S^, 51k, w^). This can be 
done with the help of the constants of motion ni, I, through equations (fZS)) . (1551) . (|571) . (pO)) . and ([57)1 . 
using the form (|72|) of the homothetic vector c. After some algebra one finds that 



Vf cx 



2\3-7u) p8(l-!ii) ^2(l+7to) 



l/(3+13«.) 



(80) 



Taking M^('^+^'^"'^/^ and replacing fl with via ([8]), and through eq. ([2g]) yields the monotone function 
Mh, which we give in the main text. Similar, but much simpler, Hamiltonian methods were used to find 
the monotone function Afcs for the non-tilted Bianchi type II models, which is also given in the main 
text. 



^This is actually the square of the monotone function in | 15| . but we find this form more convenient in the present 
context. 
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